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ABSTRACT 

We present results from the analysis of cosmic microwave background (CMB), large 
scale structure (galaxy redshift survey) and X-ray galaxy cluster (baryon fraction and 
X-ray luminosity function) data, assuming a geometrically flat cosmological model 
and allowing for tensor components and a non-negligible neutrino mass. From a com- 
bined analysis of all data, assuming three degenerate neutrinos species, we measure a 
contribution of neutrinos to the energy density of the universe, 

(68 per cent confidence limits), with zero falling on the 99 per cent confidence limit. 
This corresponds to ~ 4 per cent of the total mass density of the Universe and implies 



a species-summed neutrino mass 



0.56 



+0.30 
-0.26 



eV, 



0.2 eV per neutrino. 



We examine possible sources of systematic uncertainty in the results. Combining the 
CMB, large scale structure and cluster baryon fraction data, we measure an amplitude 
of mass fluctuations on 8ft. _1 Mpc scales of erg = 0.74^q'q7, which is consistent with 
measurements based on the X-ray luminosity function and other studies of the number 
density and evolution of galaxy clusters. This value is lower than that obtained when 
fixing a negligible neutrino mass (cr 8 = Q.86to.m)- The combination of CMB, large 
scale structure and cluster baryon fraction data also leads to remarkably tight con- 
straints on the Hubble constant, Ho — 68.4 km s _1 Mpc -1 , mean matter density, 
n ro = 0.31 ± 0.02 and physical baryon density, fl h h 2 = 0.024 ± 0.001, of the Universe. 

Key words: cosmic microwave background — cosmological parameters — dark 
matter — large-scale structure of the universe — X-rays: galaxies: clusters 



1 INTRODUCTION 

The first results from the Wilkinson Microwave Anisotropy 
Probe (WMAP; Bennett et al. 2003; Spergel et at 2003) 
have confirmed, at high statistical significance, what has be- 
come known as the 'standard' cosmological model. We ap- 
pear to live in an approximately flat, vacuum energy dom- 
inated universe, with a mean matter density Q m ~ 0.3 and 
Hubble constant Ho ~ 70 kms^ 1 Mpc~ , that was seeded by 
predominantly adiabatic, Gaussian primordial fluctuations. 

The initial papers from the WMAP team concentrated 
on results extracted from the WMAP data alone, and on 
results obtained by combining the WMAP and other cos- 
mic microwave background (CMB) data with information 
from the 2dF Galaxy Redshift Survey (2dFGRS) and obser- 
vations of the Lyman-a forest. In this paper, we present re- 
sults based on similar CMB and 2dFGRS data, but drop the 
byman-Q forest constraints. Instead, we include constraints 
provided by the baryonic mass fraction in X-ray luminous, 
dynamically relaxed galaxy clusters studied with the Chan- 
dra Observatory, and constraints obtained from the local X- 



ray luminosity function (XLF) of the most X-ray luminous 
galaxy clusters (Allen et al. 2002, 2003). 

Neutrinos are fundamental particles of the Standard 
Model of particle physics, and their masses have signifi- 
cant implications. Experiments to observe neutrino oscil- 
lations have shown there to be small, but non-zero, mass 
squared differences between three neutrino mass eigenstates 
(e.g. Fukuda et al. 1998; Ahmad et al. 2002; Eguchi et al. 
2003). If the lowest eigenstate has zero mass then these dif- 
ferences imply a mass of ~ 0.05 eV for at least one of the 
remaining two mass eigenstates. Due to the effect of the 
neutrino mass on the growth of structure in the Universe, 
cosmological observations are sensitive to the neutrino mass 
density. So far, only upper limits have been obtained from 
cosmological observations. The WMAP team (Spergel et al. 
2003) found m„ < 0.23 eV per neutrino (95 per cent confi- 
dence) from combining CMB, 2dFGRS and Lyman-Q forest 
data, assuming three neutrino species of degenerate mass 
(see Elgaroy & Lahav 2003 for a more detailed discussion of 
the assumptions). 
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Our primary motivation for this work was to compare 
and combine the constraints from the new CMB and X- 
ray cluster data, with particular attention being paid to the 
linearly-evolved amplitude of mass fluctuations. Whereas 
the CMB data probe large to intermediate length scales 
in the early universe (z ~ 1100), the XLF constraints are 
much more local (z < 0.3) and probe intermediate length 
scales. These different properties also make the combination 
of CMB and XLF data well suited to examining the contri- 
bution of neutrinos to the mass density of the universe. 

The amplitude of spatial fluctuations in the mass den- 
sity is a strong function of both scale and epoch. Con- 
ventionally, the amplitude of fluctuations is quoted on an 
intermediate length scale at the present day (z=0): ag is 
the root-mean-square matter fluctuation within 8ri _1 Mpc 
spheres, calculated using linear theory. As discussed by 
e.g. Efstathiou & Bond (1999), there are a large number 
of different (degenerate) combinations of parameters that 
can give almost identical CMB power spectra, most im- 
portantly (assuming a flat universe) the Hubble constant, 
Ho = lOO/ikms -1 Mpc -1 , the neutrino mass density, Q, v , 
and the mean matter density in units of the critical den- 
sity, fi m . Current CMB data alone cannot, therefore, pro- 
vide strong constraints on og (Spergel et al. 2003, Bridle et 
al. 2003). 

Although the normalization of the CMB power spec- 
trum on large scales is insensitive to the neutrino mass, al- 
lowing for the presence of massive neutrinos will lower the 
value of erg inferred from the CMB data. This is because 
massive neutrinos affect the growth of fluctuations on small 
scales at late times (<rg being calculated at the present day) . 
Conversely, whereas the number density of galaxy clusters 
in the local universe is sensitive to the neutrino mass, the 
value of og inferred from the XLF observations is only mildly 
affected, since the length scale probed by clusters is very 
similar to 8/jT 1 Mpc. 

Given an accurate measurement of O m , the cluster XLF 
data can yield a precise determination of erg . fi m can be ob- 
tained from measurements of the baryonic gas mass fraction, 
/g as , in large, dynamically relaxed clusters of galaxies, to- 
gether with results on the Hubble constant and mean baryon 
density of the Universe from the ratio of first to second peak 
heights in the CMB temperature angular power spectrum. 
Thus, the combination of CMB, / gas and cluster XLF data 
can together yield an estimate of the neutrino mass. 

Our approach in constraining the neutrino mass is sim- 
ilar to that described by Fukugita, Liu & Sugiyama (2000). 
A number of earlier papers have also placed limits on the 
neutrino mass using the observed slope of the galaxy (and 
by assumption matter) power spectrum obtained from the 
2dFGRS, which is sensitive to a combination of cosmolog- 
ical parameters (see Elgaroy et al. 2002, Elgaroy & Lahav 
2003, Hannestad 2003 and references therein). We show that 
including the 2dFGRS constraint is helpful in constraining 
the neutrino mass, but not essential to our analysis. 



2 METHOD 

2.1 Analysis of the CMB and large scale 
structure data 

Our analysis of CMB observations uses the WMAP temper- 
ature (TT) data for multipoles / < 900 (Hinshaw et al. 2003) 
and temperature-polarization (TE) data for I < 450 (Kogut 
et al. 2003). To extend the analysis to higher multipoles 
(smaller scales) , we also include data from the Cosmic Back- 
ground Imager (CB1; Pearson et al. 2003) and Arcminute 
Cosmology Bolometer Array Receiver (ACBAR; Kuo et al. 
2003) for I > 800. The comparison of model angular power 
spectra with the WMAP data employs the likelihood calcu- 
lation routines released by the WMAP team (Verde et al. 
2003). 

For part of the analysis, we also include large scale 
structure data from the 2dFGRS. We use the power spec- 
trum obtained from the first 147,000 redshifts of the sur- 
vey, covering scales 0.02 < fc/(/iMpc _1 ) < 0.15, where non- 
linear effects are thought to be negligible (Percival et al. 
2002). Over this range, we assume that the 2dFGRS power 
spectrum is directly proportional to the matter power spec- 
trum at z = 0, marginalizing over the constant of propor- 
tionality using a wide prior. 

Our analysis of the CMB and 2dFGRS data uses the 
CosmoMC code*. This in turn uses CAMB (Lewis, Challi- 
nor & Lasenby 2000), which is based on CMBFAST (Sel- 
jak & Zaldarriaga 1996), to generate the CMB and mat- 
ter power spectrum transfer functions, and a Metropolis- 
Hastings Markov Chain Monte Carlo (MCMC) algorithm 
to explore parameter space. We used the covariance matrix 
of the parameters to improve sampling efficiency (see Lewis 
& Bridle 2002 for more details). 

We have fitted the data using a cosmological model with 
nine free parameters: the physical dark (cold+hot) matter 
and baryon densities in units of the critical density, Qdmh 
and Qbh , the Hubble constant Ho, the neutrino mass frac- 
tion /„ = f2„/f2 dm , the recombination redshift z rcc (at which 
the reionization fraction is a half, assuming instantaneous 
reionization) , the amplitude of the scalar power spectrum 
As, the scalar and tensor spectral indexes ns and m, and 
the ratio of the normalization of tensor and scalar compo- 
nents R — At/ As, where ^4s is evaluated at k — 0.05Mpc _1 
and At is evaluated at k = 0.002Mpc _1 . We assume three 
neutrinos of degenerate (equal) mass throughout. Each pa- 
rameter had a wide uniform prior applied. We assume a flat 
geometry with a cosmological constant (f2 m + Qa=1)- 

The analysis was carried out on the Cambridge X-ray 
group Linux cluster. For our main analyses of the CMB and 
CMB+2dFGRS data, using the nine-parameter model and 
allowing for massive neutrinos and tensors, we accumulated 
at total of 10 6 correlated samples in 10 separate chains. In 
order to investigate the effects of neutrinos and tensors in 
detail, we also carried out two additional runs in which the 
neutrino density and tensor amplitude ratio were separately 
fixed to zero (/„ = and R = 0, respectively). For these 
runs, we accumulated 6 x 10 5 samples in 6 separate chains. 
We satisfied ourselves that the chains had converged by en- 
suring that consistent final results were obtained from nu- 

* http: / / cosmologist .info/ cosmomc / 
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Figure 1. The enclosed baryon fraction relative to the universal 
value as a function of radius, in units of the virial radius r v ; r , 
from the simulations of Eke et al. (1998). The simulated clusters 
have similar masses and temperatures to the clusters used here. 
The results for the most dynamically active system in the sim- 
ulations have been excluded. Beyond a radius r > 0.2r v i r , the 
simulated clusters exhibit consistent, relatively flat b profiles. At 
r = 0.25r v i r , a radius comparable to the measurement radius for 
the Chandra observations, the simulations give b = 0.824 ±0.033. 



merous small subsets of the chains. In all cases, we allowed 
a conservative burn-in period of 10 4 samples. 



2.2 Importance sampling using the X-ray cluster 
constraints 

In order to obtain improved results on cosmological pa- 
rameters we have importance sampled the CMB and 
CMB+2dFGRS MCMC results, applying further, indepen- 
dent constraints based on X-ray studies of galaxy clusters. 
Firstly, we applied constraints from the observed X-ray gas 
mass fraction, / gas , in dynamically relaxed clusters. These 
data provide one of the most simple and robust cosmolog- 
ical probes currently available, constraining a combination 
of Sim, Ho and fib- We have used the data and methods 
of Allen, Schmidt & Fabian (2002) and Allen et al. (2003). 
The data are drawn from Chandra observations of X-ray lu- 
minous, dynamically relaxed clusters spanning the redshift 
range 0.08 < z < 0.47. Consistent measurements of the mass 
profiles in the clusters are available from X-ray and gravi- 
tational lensing methods in most cases, keeping systematic 
uncertainties to a minimum. 

Each of the MCMC samples obtained from the CMB 
and CMB+2dFGRS analyses provides values for fib, fi m 
and h. Using these values, we fit the observed /gas(z) data 
with the model function (Allen et al. 2002) 



f™ od (z) - 

./gas \^) 



bfi b 



(l+0.19Vh) Sl n 



h D 



n m =i,n A 



=0/ \ 



0.5 D\ 



(1) 



where fi m = 1.0, SI a = and h — 0.5 define the reference 
cosmology used for the /gas (2) measurements. The X-ray 
gas dominates the baryonic mass content of the clusters and 
the factor (1 + 0.19aA) converts the X-ray gas mass into 



the total baryonic mass (Fukugita, Hogan & Peebles 1998). 
The quantity & is a bias factor, motivated by gasdynamical 
simulations which suggest that the baryon fraction in clus- 
ters is slightly lower than for the universe as a whole (e.g. 
Cen & Ostriker 1994; Eke, Navarro & Frenk 1998; Bialek 
et al. 2001). We use the results of Eke et al. (1998) from 
simulations of 10 X-ray luminous, massive clusters to con- 
strain b. Excluding the data for the most dynamically ac- 
tive cluster in the study of Eke et al. (1998; recall that the 
/gas data are drawn from Chandra observations of dynam- 
ically relaxed systems), the simulated clusters show consis- 
tent, relatively flat baryonic mass fraction profiles for radii 
r > 0.2r v i r (Fig 1). At r = 0.25r v i r , a radius comparable to 
the measurement radius for the Chandra observations, the 
simulations of Eke et al. (1998) give b = 0.824 ± 0.033. 

Employing a Gaussian prior on the bias factor and fit- 
ting the model described by Eq. 1, we obtain a \ 2 value for 

each of the MCMC samples. The weight of the MCMC sam- 

— 2 11 

pies is then multiplied by e x ' . Note that we assume that 
the effect of massive neutrinos on the / gas constraints is neg- 
ligible, which would be true if the neutrino to dark matter 
ratio within the galaxy clusters were typical of the universe 
as a whole. However, even in the most extreme case in which 
neutrinos completely avoid galaxy clusters, Sim in the pre- 
factor of Eq. 1 need only be replaced by Sim — Sl„. Since we 
find fi„/fi m ~ 0.04 from this work, even this extreme case 
would only increase Sim by ~ 0.01 and have negligible effect 
on our conclusions regarding the neutrino mass. 

Secondly, we have also importance sampled the MCMC 
results applying additional constraints from the observed, 
local X-ray luminosity function (XLF) of galaxy clusters. 
For this analysis, we use the results of Allen et al. (2003; 
see Fig. 3 of that paper), which are based on the XLF for 
the most X-ray luminous (Lx, 0.1-2.4 > 10 45 ergs -1 ) galaxy 
clusters within redshift z < 0.3 observed in the ROSAT All- 
Sky Survey (Ebcling et al. 2000; Bohringer et al. 2002), and 
a mass-luminosity relation determined from Chandra X-ray 
and weak lensing observations. The XLF results constrain 
the combination of Sim and as. To importance sample, we 
multiply the weight of each of the MCMC samples by the 
probability of that as, Sim parameter pair from the analysis 
of Allen et al. (2003). 



3 RESULTS 

The marginalized results on cosmological parameters ob- 
tained from the analysis of the CMB+2dF, CMB+2dF+ 
/ g as, CMB+2dF+XLF, CMB+/ gas +XLF, and CMB+2dF+ 
/gas+XLF data, using the nine parameter model with the 
neutrino mass fraction and tensor amplitude ratio included 
as free parameters, are summarized in Table 1. This table 
also lists the values of other quantities that can be derived 
from the parameter values: the matter density in units of 
the critical density, Sim {Sim = fidm + fib), the linearly- 
evolved amplitude of mass fluctuations on 8ft -1 Mpc scales, 
as, the physical neutrino density in units of the critical den- 
sity, fi„/i 2 , and the optical depth to reionization, r. 

The marginalized, posterior probability distributions 
for the parameters are also shown in Fig. 2. The results 
obtained from the additional runs with the neutrino mass 
fraction, and then the tensor amplitude ratio, fixed to zero 
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Figure 2. Marginalized probability distributions for the cosmological parameters from the analysis of the CMB+2dF (dotted) 
CMB+2dF+/ gas (dark solid) and CMB+2dF+/ gas +XLF (grey solid) data using the nine-parameter model, allowing for the presence of 
tensor components and massive neutrinos. 



are, in general, consistent with those listed in Table 1. The 
most notable exception to this is the result on as in the 
case of a negligible neutrino mass, which is discussed below. 
We note, however, that the run without tensors gives peak- 
probability values for Qbh 2 , Ho, f v , ns and ft v h 2 that are 
slightly lower (and a result on as that is slightly higher) than 
the tabulated values, although consistent within 68 per cent 
uncertainties. We also note that the limit on the ratio of the 
scalar to tensor amplitudes is constrained more tightly in 
the /„ = case (R < 0.3 at 95 per cent confidence). 



3.1 The cosmic neutrino density and amplitude of 
mass fluctuations 

Two of the more interesting results to come from the anal- 
ysis are the effects of allowing for massive neutrinos on the 
present day, linearly-evolved amplitude of mass fluctuations, 
as, and the implied contribution of such neutrinos to the en- 
ergy density of the universe, Q v h 2 . 

Fig. 3 shows the joint 68.3 and 95.4 per cent confidence 
limits on S7 m and as obtained from the CMB+2dF-r-/ gas 
data. The independent constraints from the cluster XLF 
are also shown overlaid. In the absence of massive neu- 
trinos (/„ = 0), we obtain a marginalized result on the 
amplitude of mass fluctuations from the CMB+2dF-|-/gaa 
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Table 1. The marginalized results on cosmological parameters obtained using the nine-parameter model, allowing for tensor components 
and massive neutrinos. The quoted values correspond to the peaks of the marginalized probability distributions. Error bars are 68.3 per 
cent confidence limits. Where only a limit is quoted, this corresponds to the 95.4 per cent confidence level. 



CMB+2dF CMB+2dF+/ gas CMB+2dF+XLF CMB+/ gas +XLF CMB+2dF+/ gas +XLF 



f2 dm /i 2 

a h h 2 

Ho 
fv 

A s 
n s 
JIT 


120+ 008 
U ' iZU -0.010 

0240+ 0021 
u.uz^u_ 0016 

70+ 5 

u±l 

22+1 

1 00+ 06 
1 ' uu -0.04 

0.4 ±0.6 


121+ 007 
u - lzl -0.007 

U.U2d8_ 0Q12 
68.4+ 2 -° 
0.042t° °lt 

utl 

22.0+ A 2 \ 

1 00 +oi)3 
1 - uu -0.04 

3+ - 7 


, -,.+0.008 
U - iiO -0.008 

0236+ 0018 
u.uzjd_ 0014 

71 r: + 3.5 

0.046lg : gJ 
13±g 

21. S+jg 

1 00+ 05 

4+0- 5 
U ' 4 -0.5 


n o+0. 006 
U-11B_ .006 

0236+ 0012 
u.uzod_ 0013 

67.6 +2 ' 2 

042+ 053 

20.0+fl 
1 00 +0i)3 
2 +0 ' 8 


n 1 iq+0.005 

U-ll»_ .006 
0236+ - 0010 

U.UZOD_g 0013 

68A+.H 

n^l+o- 026 
u.uoi_ 025 

12+ 4 

21.5 +2 ' 5 

0.98±g;g* 


R 


< 0.6 


< 0.5 


< 0.6 


< 0.5 


< 0.5 


n v h 2 

T 


0.28t° ; p 

0040+ ' 0035 
u.uutu_ 0035 

13+009 


0.308+° .™ 
0.74+°^ 
0045+ - 0036 

n ,0+0.06 


n 97 +0.06 

U - Zl -0.03 

72+ ' 04 
u - ,z -0.06 

0055+ 0034 
u.uuoo_ 0029 

09+ 08 


296+ 036 

69+ 04 
u - oy -0.05 

0055+ - 0070 
u.uuoo_ 0017 

0Q7 +0.15 
u - u '-0.03 


301+ 024 

70+ - 04 
— 0.04 

0059+ 0033 
u.uuoy_ 0027 

09+ 06 



0.6 



1 


1 , , , 

/ y \ 


1 




















1 


1 


1 


0.2 


0.3 
m 


0.4 



Figure 3. Joint 68.3 and 95.4 per cent confidence limits on erg 
and f2 m from the analysis of the CMB+2dF+/ gas data, allowing 
(solid curve) and excluding (dashed curve) a contribution to the 
energy density of the universe from massive neutrinos. The dotted 
curve shows the 95.4 per cent confidence limit from the analysis 
of the local cluster XLF from Allen et al. 2003 



data of erg — 0.861 o'.o? (dashed curve in Fig. 3), which is 
only marginally consistent with the cluster XLF constraints. 
When the neutrino density is included as a free parame- 
ter, however, the agreement with the XLF results is im- 
proved significantly, with a marginalized result on erg from 
the CMB+2dF+/ gas data of cr 8 = 0.74t° J 2 (solid curve in 
Fig. 3. A similar degeneracy between erg and /„ is visible in 
the analysis of pre-WMAP CMB (and other) data by Lewis 
& Bridle 2003) . We note that the effects of neutrinos on the 



cluster XLF analysis are negligible, reducing the measured 
value of as by only ~ 1.5 per cent from the no-neutrino case 
(see Section 4). 

The marginalized results on <rg from the CMB+2dF+ 
/gas data are also consistent with those obtained from other 
recent studies of the number density and evolution of galaxy 
clusters, which generally find as ~ 0.65 — 0.75 for il m = 0.3 
(e.g. Borgani et al. 2001; Reiprich & Bohringer 2002; Viana, 
Nichol & Liddle 2002; Schuecker et al. 2003; Pierpaoli et al. 
2003; Voevodkin & Vikhlinin 2003 and references therein) 
as well as some of the more recent measurements of cosmic 
shear on large angular scales (Brown et al. 2003; Jarvis et 
al. 2003 and references therein). We caution, however, that 
some other cosmic shear studies (e.g. Van Waerbeke et al. 
2002; Refregier, Rhodes & Groth 2002) have inferred higher 
values of erg ~ 0.9, for fi m ~ 0.3. 

Fig. 4(a) shows the constraints on Q, v h 2 obtained from 
various combinations of the data, as described in Table 1. 
Using the CMB±2dF±/ gas data, we measure a neutrino den- 
sity fi„h 2 = 0.0045lo;°o3i , with a 95.4 per cent confidence 
upper limit of Q. u h 2 < 0.012 (or /„ < 0.1; the same up- 
per limit is also obtained using only the CMB±2dF data.) 
For the full CMB+2dF+/ gas +XLF data set, we obtain 
ttvh 2 = 0.0059i"g; , o2 7 , with the region between Q. v h 2 = 
and the corresponding equi- probability upper value for Q, v h 2 
containing 99 per cent of the probability. Finally, for the 
CMB±/ gas +XLF data, we find VL v h 2 = 0.0055lg;g o i?, with 
£l u h 2 > at more than 99 per cent confidence. We thus see 
that including the 2dFGRS constraint is helpful in the anal- 
ysis, but not essential to our detection of a non-zero neutrino 
mass. Our most probable result on the neutrino density from 
the CMB+/ gas +XLF and CMB±2dF+/ gas +XLF data sets 
corresponds to ~ 1.3 per cent of the critical value required 
for closure. 

Fig. 4(b) shows how the results on Q. v h 2 are modified 
when we exclude the presence of tensor components. In this 
case, the full CMB±2dF+/ gas +XLF data set gives Q, v h 2 = 
0.0034+Q o^g. Although the most probable value for n„h 2 
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Figure 4. (a) Marginalized probability distribution for Q, v h? from the analysis of the CMB+2dF+/ gas (dotted curve), CMB+/ gas +XLF 
data (grey curve) and CMB+2dF+/ gas +XLF (dark, solid curve) data. In all cases we have allowed for the presence of tensor components 
(b) The dashed curve shows the results from the CMB+2dF+/ gas +XLF data from the analysis without tensor components. The dark, 
solid curve is the same as in (a). 



is reduced when we exclude tensor components, we still find 
Q v h 2 > at 93 per cent confidence. 

The results on Q, u h 2 can be converted into a con- 
straint on the neutrino mass, summed over species, 
using the relation Q. u h 2 — ^.mj/94eV. For the 
full CMB+2dF+/ gas +XLF data set, we find J2i m » = 
0.56t o j°eV. 

Our result on the contribution of neutrinos to the en- 
ergy density of the universe differs slightly from that re- 
ported by the WMAP team, Q. u h 2 < 0.0076 (95 per cent 
confidence limit; Spergel et al. 2003). The primary reason 
for this difference is the inclusion of Lyman-a forest data 
in the WMAP team's analysis (see discussions in Elgaroy 
& Lahav 2003 and Hannestad 2003) and the inclusion of 
cluster XLF data in ours. Our results are, however, consis- 
tent with those of 2dFGRS team: J2i m * < 2 - 2 eV ( E lg ar °y 
et al. 2002; Hannestad 2003) as well as those of Lewis & 
Bridle (2003) and Slosar et al. (2003) from the analysis of 
pre- WMAP CMB, 2dFGRS and other data using MCMC 
techniques. 



3.2 The Hubble constant and mean matter 
density 

Our results on the Hubble constant are in good agreement 
with, but more tightly constrained than, that from the Hub- 
ble Key Project (Ho = 72 ± 8kms _1 Mpc" 1 ; Freedman 
et al. 2001). Including the neutrino energy density as a 
free parameter, our analysis of the CMB+2dF data gives 
Ho = 70±5 km s" 1 Mpc" 1 . The introduction of the Chan- 
dra / gas data improves this constraint significantly, to Ho = 
68.412-Okms- 1 Mpc -1 . For the full CMB+2dF+/ gas +XLF 
data set, we find H = 68.4t{\l kms" 1 Mpc" 1 . 

The introduction of the / gas data also leads to signifi- 



cant improvements in our measurements of the mean mat- 
ter and baryon densities. As can be seen from Table 1, the 
combination of CMB+2dF data gives a ~ 20 per cent uncer- 
tainty on fi m and an eight per cent uncertainty on Q,\>h 2 . The 
introduction of the / gas data improves these constraints to 
seven per cent on fl m and five per cent on Q^h 2 . The further 
inclusion of the XLF data makes little difference to the re- 
sults, providing us with final answers of Q m — 0.301 ±0.024 
and fl^h 2 = 0.0236 ± 0.0012. The improvement in the con- 
straints on Ho and fi m obtained with the introduction of the 
Chandra / gas data can be clearly seen in Fig. 5. 

We note that the tight constraints on Ho, Q m and 
flbh 2 are primarily due to the combination of CMB and 
/ gas data. From an analysis of these two data sets alone, we 
find H = 66.4±^ kms" 1 Mpc" 1 , Q m = 0.3041° ^ and 
Q b h 2 = 0.0234j£°jg£. 



4 DISCUSSION 

The main result of this paper is the tentative detection of a 
non-negligible neutrino mass. Recent years have seen signif- 
icant progress in our understanding of neutrinos from par- 
ticle physics. The Large Electron Positron Collider (LEP) 
and Stanford Linear Accelerator Centre Linear Electron- 
Positron Collider (SLC) showed there to be three light, ac- 
tive neutrino species (electron, muon, tau). Solar and at- 
mospheric neutrino experiments such as Super-Kamiokande, 
the Sudbury Neutrino Observatory and Kamland (Fukuda 
et al. 1998; Ahmad et al. 2002; Eguchi et al. 2003) have since 
provided strong evidence for oscillations between these neu- 
trino species, and determined two of the mass-squared dif- 
ferences to be Am s 2 ol ~ 7 x 10" 5 eV 2 and Am atm ~3x 10" 3 
eV 2 . Given the relatively small mass-squared differences, the 
result on the species-summed neutrino mass presented here, 
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£\ rm = 0.56±° |g eV (for the CMB+2dF+/ gas +XLF data) 
implies an approximately degenerate neutrino mass. Writ- 
ing rm = Nm v and setting N = 3 (assuming an ab- 
sence of further, heavy sterile neutrino species), we obtain 
m„ = 0.19tg'Jg eV. Note, however, that our results are also 
marginally consistent with indications from the Liquid Scin- 
tillator Neutrino Detector experiment (Aguilar et al. 2001) 
for the presence of sterile neutrinos. 

Our results are consistent with the upper limit on the 
electron neutrino mass, m v < 2.2 eV (95 per cent confi- 
dence limit), from laboratory tritium beta decay end-point 
experiments (Weinheimer et al. 1999; Bonn et al. 2002). 
We note that our results are also consistent with the neu- 
trino mass inferred from Z-burst models (Fargion, Mele & 
Salis 1999; Weiler 1999) for ultrahigh energy cosmic rays, 
m v = 0.26±o 1° eV (Fodor, Katz & Ringwald 2001; assum- 
ing an extragalactic origin for the cosmic rays), some in- 
dications from neutrinoless double beta decay experiments, 
0.11 eV < m v < 0.56 eV (Klapdor-Kleingrothaus et al. 2002; 
although see also Aalseth et al. 2002; Feruglio, Strumia & 
Vissani 2002) and the neutrino mass range inferred from 
anthropic arguments by Tegmark & Vilenkin (2003). 

As was discussed in Section 1, our result on the neutrino 
mass density originates primarily from the fact that the XLF 
provides a robust constraint on as for a given value of f2 m , 
while the CMB data predict as as a function of the neutrino 
mass. The primary effect of the introduction of massive neu- 
trinos is to suppress the formation and growth of structure 
on small scales (k < 0.026(m„/l eV)°' 5 fi^ 5 h Mpc" 1 ; Hu, 
Eisenstein & Tegmark 1998). The large momenta of such 
neutrinos prevents them from clustering with the cold mass 
components and they stream freely. The suppression of the 
matter power spectrum as a function of scale in the case of 
massive neutrinos, relative to the massless neutrino case, is 
shown in Fig. 6. 

Since the CMB data normalize the power spectrum pri- 
marily on large scales, the amplitude of fluctuations on (rel- 
atively small) Sh^ 1 Mpc scales, as, inferred from the CMB 
data will always be smaller in the presence of massive neutri- 
nos (~ 15 per cent smaller for Q. v h 2 — 0.0059; Fig 6). In con- 
trast, the value of as inferred from the cluster XLF analysis 
is relatively insensitive to massive neutrinos, since clusters 
measure the power on scales similar to 8 h~ Mpc. (The mean 
length scale corresponding to the virial masses of the clus- 
ters in the XLF study of Allen et al. 2003 is 13.6 /i _1 Mpc.) 
For f]„/i 2 = 0.0059, the value of as inferred from the XLF 
study is reduced by only ~ 1.3 per cent when accounting for 
the effects of the neutrinos. 

We note that on scales significantly smaller than those 
of clusters, accounting for the presence of neutrinos at 
the measured level will increase the inferred value of as 
by a small amount (~ 2.1 per cent for measurements on 
2/i~ 1 Mpc scales, for Q v h 2 = 0.0059). This may be relevant 
to studies of cosmic shear and the Lyman— a forest. Note, 
however, that non-linear corrections to the matter power 
spectrum also become important on such scales {e.g. Smith 
et al. 2003). 

In what follows, we examine the most significant as- 
sumptions and systematic uncertainties affecting our result 
on the neutrino mass. These primarily concern the data sets 
used and the parameters varied in the analysis. 




Figure 5. The joint 68.3 and 95.4 per cent confidence limits on 
Ho and f2 m from the analysis of the CMB+2dF (dotted curve) 
and CMB+2dF+/ gas (solid curve) data. We have allowed for the 
presence of tensor components and massive neutrinos in the anal- 
ysis. 
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Figure 6. The suppression of the matter power spectrum as 
a function of scale due to massive neutrinos, relative to the 
massless neutrino case. The solid curve shows the suppres- 
sion for f2„/i 2 = 0.0059, the peak-probability value from the 
CMB+2dF+/ gas +XLF analysis, allowing for tensor components. 
The dashed curve shows the results for Q„h 2 = 0.0045, the result 
from the CMB+2dF+/ gas data. All other parameters were fixed 
to their best-fit values (Table 1). 
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• Reliability of XLF constraints. As was mentioned 
above, the result on as as a function of $l m from the clus- 
ter XLF study is crucial to our measurement of the neu- 
trino mass. Systematic uncertainties in the XLF analysis 
are discussed by Allen et al. (2003). The main uncertainty 
in that work lies in the mass-luminosity relation used to con- 
vert model mass functions into X-ray luminosity functions, 
which are then compared with the observations. Allen et al. 
(2003) find a» = 0.69 ± 0.04 for fi m = 0.3, accounting for 
uncertainties in the normalization, slope and scatter of the 
mass-luminosity relation. Replacing the Allen et al. (2003) 
XLF constraints with a Gaussian probability density func- 
tion (pdf) with erg = 0.69 ± 0.04 leads to similar results on 
the neutrino mass. 

We have examined the effects on the measured neu- 
trino mass of introducing additional, systematic offsets 
in the mass-luminosity relation. For example, replacing 
the Allen et al. (2003) XLF constraints with a Gaus- 
sian pdf with as = 0.77 ± 0.04 (which corresponds to 
raising the normalization of the mass-luminosity curve 
by a factor ~ 1.3) reduces the neutrino mass inferred 
from the CMB+2dF+/ gas + (replaced) XLF data to m v = 
0.14+Q'Jg eV. Conversely, lowering the normalization of the 
mass-luminosity relation by a similar factor using a Gaus- 
sian pdf with erg = 0.62 ± 0.04 increases the neutrino 
mass to m v = 0.26lg'°g eV. (Note that the inclusion of the 
CMB+/gas data maintains fi m ~ 0.3.) 

The results on the neutrino mass are not sensitive to 
moderate changes in the uncertainty in as from the XLF 
data. Approximately doubling the uncertainty in as using 
a Gaussian pdf with as = 0.69 ± 0.08 still gives m„ = 
0.19 ± 0.10 eV with Q. v h 2 > at more than 99 per cent 
confidence for the CMB+2dF+/ gas + (replaced) XLF data 
(allowing for tensor components). Only if the uncertainty 
in as is quadrupled from its nominal value of 0.04 does the 
detection of a non-zero neutrino mass drop below the 95 per 
cent confidence level. 

• Reliability of CMB polarization data. The CMB 

temperature-polarization correlation power spectrum helps 
to constrain the optical depth to reionization, r. The most 
probable values for r found from this analysis and by the 
WMAP team (Kogut et al. 2003; Spergel et al. 2003) are 
higher than were expected prior to the release of the WMAP 
results. If we were to exclude the CMB polarization data, 
then the CMB would constrain only the degenerate combi- 
nation (7ge _T . Since this product decreases with increasing r, 
the CMB+2dF+/ gas constraint could also be matched with 
the XLF constraint on as if the value of r were lowered. For 
example, if we constrain 0.04 < r < 0.07 in the analysis of 
the CMB+2dF+/ gas data, then we obtain a» = 0.73 ± 0.07. 
In this case, when combining with the XLF data, a negligi- 
ble neutrino mass falls on the ~ 96 instead of 99 per cent 
confidence limit. 

We note that the measured optical depth to reioniza- 
tion, r — 0.09^ 05 from our default analysis of the full 
CMB+2dF+/ gas +XLF data set, is lower, and has larger un- 
certainties, than the value quoted in the abstracts of the 
WMAP papers (0.17 ± 0.04; Kogut et al. 2003; Spergel et 
al. 2003). Taken together with the constraints on ns, our re- 
sults may be slightly easier to explain within the context of 
standard UV reionization models, even allowing for the sup- 



pression of power on small scales due to the free-streaming 
action of massive neutrinos. 

• Reliability of the 2dFGRS. The 2dFGRS constraints 
used in this paper rely on the assumption of linear bias. This 
issue was discussed by Elgaroy & Lahav (2003). However, 
the main point here is that our result on the neutrino mass 
is largely unchanged when removing the 2dFGRS constraint 
(Section 3.1). 

• The effect of tensors. Including the presence of pri- 
mordial gravity waves in our analysis increases power on 
the largest scales. This is reflected by the fact that a fit 
to the CMB+2dF+/ gas -|-XLF data in the presence of ten- 
sor components gives a slightly higher scalar spectral index, 
ns — 0.981q;q3, than is obtained in the absence of tensors 
ns = 0.95 ± 0.02. The higher scalar spectral index in the 
presence of tensors slightly increases the power on interme- 
diate scales, relative to the largest scales and, therefore, also 
slightly increases the neutrino mass relative to the no-tensors 
case. However, the detection of a non-zero neutrino mass in 
the absence of tensors remains significant (Section 3.1). 

We note that the most probable scalar spectral index 
obtained from our analysis is consistent with standard in- 
flation models. The ratio of the tensor to scalar compo- 
nents is also consistent with single-field, slow-roll inflation 
(R = -8n T ). 

• The effect of a running spectral index. There has 
been recent speculation about a non-power law primordial 
power spectrum. Evidence for a running spectral index was 
found by Spergel et al. (2003), when the WMAP and other 
CMB data were combined with 2dFGRS and Lyman-a for- 
est constraints. It is interesting to consider whether our con- 
straint on a significant neutrino mass could be negated by al- 
lowing for a non-zero value of n run . For the CMB+2dF+/ gas 
data, we find an n run value consistent with that of Spergel et 
al. (2003), but with wider error bars (assuming a negligible 
neutrino mass). This is mainly due to the fact that we do 
not use the Lyman-a forest data in our analysis. However, 
we find that although the constraint on as is widened when 
including n run as a free parameter, the lower limit on as 
is hardly changed. This acts in an opposite sense to allow- 
ing freedom in the neutrino mass and does not improve the 
agreement between CMB+2dF+/ gas and XLF constraints 
shown in Fig. 3. 

• The effect of quintessence. We can lower the value of 
as allowed by the CMB data if we increase the dark en- 
ergy equation of state parameter, w, above the cosmological 
constant-like value of w = — 1 (see e.g. Bridle et al. 2003b). 
Repeating the analysis of CMB data, but replacing the free 
neutrino mass by a free equation of state parameter, we ob- 
tain w ~ —0.6 for as ~ 0.7. However, such a solution is 
disfavoured by the combination of CMB+/ gas data, which 
prefer a value for w ~ — 1. Our result on quintessence is in 
line with supernovae studies, which also require w ~ — 1 for 
n m ~ 0.3 (Tonry et al. 2003). 

In future work, we will examine in more detail the con- 
straints that the combination of CMB and X-ray cluster 
data can place on the dark energy density and w. 
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